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Abstract 

Including the vacuum effects, the compressional properties of nuclear matter 
are studied in the cutoff field theory. Under the Hartree approximation, the 
low-energy effective Lagrangian is derived in the framework of the renormal- 
ization group methods. The coefficients are determined in a way where the 
physical results hardly depend on the value of the cutoff which is conveniently 
introduced into the theory. It is shown that, to reproduce the empirical data 
of the nucleus incompressibility, the compressibility of the nuclear matter is 
favorable to be 250~350MeV. 
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1. Introduction 



In recent two decades, nuclear matter has been studied in the frame- 
work of quantum hadrodynamics (QHD). The meson mean-field theory for 
nuclear matter [1] has made successful results to account for the saturation 
properties at the normal nuclear density. Following to those successes, many 
studies and modifications are done in the relativistic nuclear models. One 
of those modifications is inclusion of vacuum fluctuation effects, which cause 
divergences in physical quantities as they are naively calculated. Chin [2] 
estimated the vacuum fluctuation effects in the Hartree approximation by 
using the renormalization procedures, and found that the vacuum fluctua- 
tion effects make the incompressibility of nuclear matter smaller and closer 
to the empirical value than in the original Walecka model. However, it be- 
comes more difficult to do the renormalization as the model becomes more 
complicated, since the renormalization procedures need analytical studies to 
some extent. Numerical studies using the cutoff field theory may be useful. 

Furthermore, the relation between QHD and the underlying fundamental 
theory, i.e., QCD, is an open question. One may wonder whether QHD is 
valid in very high-energy scale or not. If QHD is valid only under some energy 
scale, it is natural to introduce a cutoff or a form factor into the theory. One 
may introduce the cutoff [3] or the form factor [4] to avoid the instability 
of the meson propagators in the random phase approximation (RPA) [5]. 
Cohen [6] introduced the four dimensional cutoff into the relativistic Hartree 
calculation and found that the vacuum energy contribution may be somewhat 
different from the one in the renormalization procedures, if the cutoff is not 
so large. 

In the recent paper [7], we have studied the nuclear matter properties in 



2 



details using the cutoff field theory of Cohen and show that the properties 
of nuclear matter may be somewhat different from those described in the 
ordinary renormalization procedures if the cutoff is not so large. In sec. 3 of 
the same paper, the cutoff field theory of the nuclear matter is reformulated 
in the framework of the renormalization group methods [8] [9] [10] and the fifth 
and sixth order terms of a-meson self-interactions in the effective Lagrangian 
are shown to be important in the cutoff field theory. We have also studied the 
meson masses at finite baryon density [11] and the vertex corrections [12] in 
the framework of the cutoff field theory and renormalization group methods. 
In this paper, we studies the compressional properties of the nuclear matter 
in detail by using the method presented in sec. 3 of [7]. 

This paper is organized as follows. In section 2, we review our method 
to study the vacuum effects in the cutoff field theory. In section 3, we show 
the numerical analyses of the compressional properties of the nuclear matter. 
The section 4 is devoted to the summary. 

2. Cutoff Field Theory for Nuclear Matter 

At first, we start with the following renormalizable Lagrangian of a-u 
model together with a regulator that truncates the theory's state space at 
some large A. 

L = — M + g s (f) — g vli y^ + \d^d^ - l -m\<f - U{<j>) 

W) = E cmt, (i) 

n=0 
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where ip, 0, V^, M, m s , m v , g s , and g v are nucleoli field, cr-meson field, uo- 
meson field, nucleon mass, a-meson mass, u;-meson mass, a-nucleon coupling, 
and u;-nucleon coupling, respectively. The C n are constant parameters which 
are adjusted to reproduce the physical conditions as explained below. The 
Lagrangian (1) is valid only in the region of the energy scale which is smaller 
than A. 

In the relativistic Hartree approximation with the cutoff A, the one-loop 
contribution to the a effective potential is given by [6] [7] 

IWM*,A) = -/ ^^r(log(^±^))0(A 2 -A;|), (2) 

where M* = M — $ = M — g s < >, is an arbitrary scale parameter with 
dimensions of mass, G is the step function, and the subscript E denotes that 
the momentum with it is written in the Euclidian notation. (We remark that 
/^-dependence appears only in the density-independent term in the effective 
potential (2). So the /i-dependence of the energy density e also appears 
only in the density-independent term of e and disappears when we choose 
the coefficient C to ensure e = at zero density, as is described below. ) 
We choose parameters C n to reproduce the following conditions as in the 
ordinary renormalization procedures [2] [13]. 

-^[Um + U^ loop m |* =0 = 0. (n = 0, 1, 2, 3, 4) (3) 

Note that the different conditions give the different physical results as is 
pointed out by Heide and Rudaz [14]. 
The conditions give 

1 1 d n 

C n = --Ut) oop (0) = --—UW^l**, (n = 0, 1,2,3,4). (4) 

The condition (4) with n = ensures that energy density e of the system 
becomes zero at zero baryon density (p = ) and remove the /z-dependence 



4 



of e. The condition (4) with n — 1 ensures that the scalar density p s of the 
nucleons becomes zero at p = and M* = M at p = 0. The condition (4) 
with n = 2 ensures that the physical mass of a- meson is m s . The conditions 

(4) with n = 3 and n = 4 mean that the effective cubic and quartic self- 
couplings of cr-meson vanish. 

After the C n are determined, the scalar part of the nucleon self-energy is 
calculated by the equation of motion (EOM) for cr-meson, i.e., 

3(1 = 0. (5) 

The vector part of the nucleon self-energy is calculated by the equation of 
motion of cj-meson as usual [1] [2] [13]. If we consider the limit A — > oo, 
the cutoff model is equivalent to the ordinary model in the renormalization 
procedures [2] [13]. 

In the Hartree approximation, the total energy density of the system is 
given by the following equation. 

e(i „ A) = ^ f dk z \/WTW> + ^ + 

+[/(*) + U l - loor (M', A), (6) 

where C s = g s M/m s , C v = g v M/m v , g is the degeneracy factor (g = 4 in the 
symmetric nuclear matter ) and hp is the Fermi momentum. Note that eq. 

(5) is equivalent to 

&r(fc F ,M-$,A) de(kp, M*, A) = 

9$ dM* ' 1 ' 

In the sec. 2 of ref. [7], using the eqs. (6) and (7) with the coefficients (4) 
and the given value of the cutoff A, we have done self-consistent calculation 
in the Hartree approximation. However, here we use the method which is 
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presented in the sec. 3 of ref. [7] and review the derivation of the low-energy 
effective Lagrangian in the framework of the renormalization group method 
[8] [9] [10]. 

If we do not know the value of A, the cutoff which we introduce into the 
theory (below we call this cutoff A' ) may be different from the true cutoff A 
which is the limiting energy scale of the theory. Suppose A' is smaller than 
A. (Of course, A' should be larger than the energy scale of the physics in 
which we are interested. ) 

In the renormalization group methods [8] [9] [10], we require that the phys- 
ical quantities do not depend on A', although the physical quantities depend 
on A. 

To achieve this, we estimate the contributions which is needlessly dis- 
carded by introducing the cutoff A' which is smaller than A. In the case of 
the vacuum energy potential, it is given by 

AU = -AL = Ui-i oop (M*, A 2 > k% > A' 2 ) 

= U^ loop (M*, A) - U^ loop (M*, A') 

9 f A „ f kl + M* 2 ^ gA' 4 /-a/A' ,x 2 + y\ , n . 

where x = fog/A' and y = M* / A' . If M* / A' < 1, we could expand the last 
line of eq. (8) around y 2 = 0. 

i er(AU) , ^ l d*"(AU) M-* 2m 

Therefore, if we use the new cutoff A', we must add AL to the Lagrangian 
as the effective potential, to keep that physical qunatities do not depend on 
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A'. L + AL is the low energy effective Lagrangian under the energy scale 
A' in the Hartree approximation. With this effective Lagrangian, the energy 
density is obtained by 

e(k F , AT, A') — AL = e(k F , M*, A') + AU, (10) 

which is equivalent to e(k F , M* , A). It should be remarked that, since the 
terms $ n with n < 4 can be absorbed in U(<&) which is phenomenologically 
determined, only the terms with n > 5 in eq. (9) are essentially new. 

Since the coefficients of the expansion (9) are order (A') 4 , the m-th term 
of (9) has the (1 / A') 2 " 1-4 order contribution if we treat A as the same order as 
the A'. In actual calculations, we truncate the Taylor expansion (9) at some 
finite maximum m. In that case, the calculated results include 0((1/A') 2m ~ 2 ) 
errors. To get higher accuracy, higher order terms in the Taylor expansion 
(9) are needed. 

In the ordinary renormalization procedure, only the terms $ n (n < 4) 
should be determined phenomenologically, since the limit A' — > oo is taken. 
However, if A is finite and A' can not got to infinity, the errors of the or- 
der (1/A') 2 remain. According to Lepage's proposal in QED case [10], we 
proposed that not only the term $ n in < 4) but also the higher terms of $ 
should be determined phenomenologically to remove the errors which arise 
from the finiteness of A'. For an example, we should determine the coeffi- 
cients for the term (/ < 6), if we want the results with 0((1/A') 4 ) error. 
This is equivalent to use 

V = L - C 5 $ 5 - C 6 $ 6 (11) 

as the effective Lagrangian with the cutoff A', instead of the original La- 
grangian (1). The coefficient C5 and C% are determined phenomenologically 
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as well as C n (n < 4). In that case, the total energy density of the system is 
given by 

+[/($) + C/i_ ioop (M*, A') + C 5 $ 5 + C 6 $ 6 . (12) 

The error and the A'-dependence of this energy density are the order of 
(1/A') 4 . Similarly, by adding more higher terms of to the right-hand side 
of eq. (11), we can get the results, the error and the A'-dependence of which 
are the order of (1/A') 2 - 7-2 , where j = l ma x/2- ( Note that l max should be a 
even number. ) 

What is the merits of using the new Lagrangian V with A' instead of 
the original Lagrangian L with the original cutoff A? The main merits are 
following. 

(1) By using the new extended Lagrangian L f , the contributions beyond the 
approximations used in the calculations are also reincorporated phenomeno- 
logically into the calculations as well as the contributions which are discarded 
in using the convenient cutoff A' instead of A. 

(2) In the ordinary cutoff theory as is described in section 2 of ref. [7], the 
results may depend on the details of the regulator which is used in the calcu- 
lations, even if the value of A are determined phenomenologically. However, 
in the new method, such a dependence on the details of the regulator is 
removed order by order, by determining the coefficients of the Lagrangian 
phenomenologically, as well as the dependence on A' itself. 

(3) In the Hartree approximations, the momentum dependence of the physical 
quantities is not calculated. However, e.g., in the calculations of the vertex 
corrections [12], by introducing new cutoff A', it can be omitted consistently 
to calculate the momentum dependence of the quantities in high momentum 
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region {k 2 E > A' 2 ) in which we are not interested. This makes the calculations 
much easier. 

It should be also emphasized that the effective Lagrangian (11) has the 
" nonrenormalizable" terms, although the original Lagrangian is "renormaliz- 
able" . Of course, this causes no difficulty since we use (11) with the cutoff A'. 
On the other hand, if the "nonrenormalizable" Lagrangian (11) is the original 
Lagrangian at A, we must determine C 5 and C 6 phenomeno logically as well as 
g s , g v , C , Ci, C 2 , C 3 and C 4 . This indicates that, in the cutoff field theory, 
it is difficult to know whether the original Lagrangian is " renormalizable" 
or "not", unless A is exactly known. In this context, "nonrenormalizable" 
Lagrangian is natural and useful in the cutoff field theory. 

One may wonder whether the expansion (9) well converges if A and A' are 
not so large as is expected in the nuclear physics. In fact, the direct calcula- 
tion of the expansion (9) shows that the expansion does not well converges 
if A' < 5GeV [7]. However, we can re-expand (9) as 

Comparing (9) with (13), it is easily seen that Z-th term in (13) is order of 
(1/A') 2j,_4 , where j is the integer part of (I + l)/2. However, different from 
the expansion in (9), besides the contribution of the order of (1/A') 2j, ~ 4 , all 
higher order contributions which are proportional to are also included. 
It has been shown that the expansion (13) well converges even if A' is not 
so large [7]. It should be remarked that the expansion (13) is automatically 
used if the coefficients of AU are determined phenomenologically. Therefore, 
we can use the expansion (13) in nuclear physics. 

We also remark that the coefficients such like C5 and Cq are not needed 
to be observed directly. It can be also determined by the physical quantities 
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at the finite density as is in the cases of the other phenomenological models 
(See, e.g., [15]). For an example, if the effective nucleon mass M* and the 
incompressibility K are given in addition to the saturation conditions, C 5 
and Cq can be determined. Some examples of the parameter sets are already 
shown in ref. [7]. In section 3, we analyze the compressional properties of 
the nuclear matter in detail using the effective Lagrangian (11). 

3. Compressional Properties of Nuclear Matter 
in Cutoff Field Theory 

In this section, we study the compressional properties of the nuclear mat- 
ter in detail using the cutoff field theory described in the previous section. 
There are eight parameters in our model. If we require the conditions (3) as 
in the ordinary renormalization procedure, four parameters, i.e., C s , C v , C 5 
and Cq remain as the free parameters of the model. If the effective nucleon 
mass is given at the normal density p , C v is determined by the relation [16] 



which is derived from the Hugenholtz-van Hove theorem [17] at the normal 
density po- The subscript denotes that the corresponding physical quantity 
is the one at the normal density p Q . By putting C 2 U = in eq. (14), it is 
also seen that the upper limit of the effective nucleon mass at the normal 
density is 0.944M. 

The other three parameters are determined by giving the saturation con- 
dition and the value of the incompressibility 




(14) 



K = 9 Pl 



2 d 2 E l 



(15) 



dp' 



2 \P=Poi 
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where Ej, = e / p — M. 

Below, we use E b0 =-15.75MeV as the binding energy at p = 0.15fm -3 . 
We treat Mq and K as variable inputs. 

Using these inputs, we calculate the skewness coefficient which is defined 

by 

In our definition of K', large K' means that the equations of the state becomes 
hard at high densities. 

In fig. 1, we show the A' dependence of K' . In these calculations, first we 
search the parameter sets of C s , C v , C 5 and C 6 to reproduce the saturation 
conditions and the given Mq and K at A' =1.5GeV. After that, fixing C s 
and C v , we vary A' and determine C 5 and C 6 to reproduce the saturation 
conditions. Using new A' with new C 5 and C 6 , K' is calculated. It is expected 
that the effects of changing A' are almost canceled by the effects of the change 
of C 5 and C 6 , if the effective Lagrangian (11) is valid. 

As is seen from these figures, the A' dependence of K' is very small when 
M* = 0.75M and 0.9M. In the case of M* = 0.6M, there is the weak A' 
dependence. This indicates that the higher terms of $ in the AL, which 
are neglected, may become important when $ is large, i.e., Mq is small. 
However, even in these cases above, the dependence is only the order of 10 
percents of K' itself. The A' dependence is well removed by the corrections 
of the order (1/A') 2 in (11). Therefore, we fix A' = 1.5GeV below. 



Fig. 1 (a), (b) and (c) 



In fig. 2, we show the K — K' relation with fixed Mq. From the figure, 
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it is seen that, in the case of Mq < 0.8M, K' becomes larger as K becomes 
larger. In the case of Mq = 0.9M, K' decreases in the region 180MeV <K. 
This decrease of K' is related to the rapid decrease of C 5 . The C 5 becomes 
negative in the case of Mq = 0.9M and K> 250MeV, while in the other cases 
C 5 is positive. 



Fig. 2 



The K' is related to the coefficients of the leptodermous expansion [18] 
[19] of nucleus incompressibility K{A, Z) as follows. 

K(A,Z) = K+K sf A- 1/3 + K vs I 2 + K c Z 2 A~ 4/3 + --- ; I = 1-2Z/A, (17) 

where the coefficients K s f, K vs and K c are surface term coefficient, volume- 
symmetry coefficient and Coulomb coefficient, respectively. In the scaling 
model [18] [19], K c is related with K and K' via 

K <=-rt-K + * (18) 

where e is the electric charge of proton and ^f r o = VPo- Pearson [19] 
pointed out that there is strong correlation between K and K c which are de- 
termined from the experimental information of the giant isoscalar monopole 
resonances, although, at present, K and K c are not determined uniquely. ( 
He pointed out that K = 120 ~ 351MeV is available to fit the data. ) The 
relation which Pearson found is shown in table I [19]. The K vs which has 
also strong correlation with K is shown in the table, too. The similar results 
are also gotten by Shlomo and Youngblood [20] . 
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Table I 



Putting the value of K' into eq. (18), we calculate K c as the function of 
K. The results are shown in fig. 3. In the case of Mq < 0.8M, K c becomes 
smaller or more negative as K becomes larger. In these cases, the calculated 
values are close to the empirical one at each K. In the case of Mq = 0.9M, K c 
increases in the region K > 220MeV. This increase is related to the decrease 
of K' in fig. 2, however, the value of K at the starting point of this increase 
of K c becomes somewhat larger than the one at the starting point of the 
decrease of K' by the factor K in the denominator of the right-hand-side of 
eq. (18). 



Fig. 3 



Next we search the parameter sets which reproduce the empirical value 
of K and K c in table I. If we require Cq > 0, which ensures the stability of 
the solution of EOM of <r-meson, the allowed region of K is restricted. The 
results are summarized in table 11(a). If this condition are required, only 
K = 250 ± 50MeV is allowed. In all allowed cases, M * > 0.85M. 



Table 11(a) and (b) 



If the negative Cq is not forbidden, some other parameters sets are al- 
lowed. Such parameter sets which reproduce the mean value of K c in table 
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I are summarized in table 11(b). In contrast to the cases with Cq > 0, 
M* < 0.85M in table 11(b). Naturally, the solution of EOM of a-meson 
becomes unstable at the limit $ — > oo in the case of negative C 6 . However, 
since we have neglected the higher terms of $ in AL, the large <3> region may 
be out of scope of this approximation. In fig. 4, we show the derivative e' 
of total energy density with respect to $. The zero point of Jj| in the figure 
is the solution of eq. (7). The solution is stable if e' is negative for $ < $ s 
and positive for $ > $ s , where $ s is the solution. From the figure, it is 
seen that the solutions of the parameter sets PS6, PS7 and PS8 are very 
unstable. However, the solutions of the parameter sets PS9 and PS10 are 
stable in the region $ < M. Therefore, below, we keep the PS9 and PS10 as 
the candidates of the realistic parameter sets. In these cases, K lies in the 
region 250~300MeV. 



Fig. 4 



In the scaling model, volume-symmetry coefficient in eq. (17) is also 
related to K and K', i.e., 

K vs = K sym -L(9^ + Q), (19) 



where 



°~dp =p0 ' s y m = ^°~^2"' P=P0 ' a ° 4 = 2^9/? 3=0 ' 

In eq. (20), P3 — Pp — Pn, and p p and p n are proton density and neutron 
density, respectively. It is well known that p-meson contribution is impor- 
tant in the calculations of these quantities which characterize the symmetric 
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properties of the nuclear matter [21] [13]. The p-meson contribution to the 
energy density in the Hartree approximation is given by [21] [13] 

" = ^-^* = ^A (21) 

where m p and g p are p-meson mass and p-nucleon coupling, respectively, and 
C p = g p M/m p . We choose C p to reproduce the empirical value of symmetry 
energy, = 30.0MeV and calculate L, K sym and K vs including the p-meson 
contribution (21). The results are also summarized in table II and fig. 5. 
The empirical value of sets of (K, K vs ) [19] [20] are also shown in fig. 5. From 
the figure, we see that only parameter set PS4 can well reproduce Pearson 
data of K, K c and K vs at the same time. Among the other parameters set, 
PS10 gives K vs which is close to the empirical value. 



Fig. 5 



In table 11(b), for negative C 6 we only show the parameter sets which 
reproduce K and the corresponding mean value of K c in table I. If we use the 
upper bound of K c , we get other parameter sets. Among them, C 2 = 358.11, 
C 2 = 270.07, C 5 M = 0.0025532, C 6 M 2 = -0.0018579 and C 2 p = 66.135 
which are chosen to reproduce K = 350.0MeV, K c = -7.274 + 2.06MeV and 
a 4 = 30.0MeV, give M* = 0.5442M and K vs = -357.8. The value of K vs has 
good agreement with the empirical data in table I. We call this parameter set 
PS11. This parameter set has the stable solution of EOM of a-meson in the 
region $ < M. On the other hand, the parameter sets which reproduce the 
lower value of K c in the table I has a unstable solution of EOM of cr-meson. 
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Finally, we show the p-dependence of the binding energy and the effective 
nucleon mass by using the parameter sets PS1~5, 9, 10 and 11 , which have 
the stable solution of EOM of a-meson in the region $ < M. The results 
are shown in fig. 6 and fig. 7. From fig. 6, it is seen that the equations of 
state is much more harder at high densities in the cases of PS9, PS10 and 
PS11 than the other cases with C 6 > 0. Samely from fig. 7, we see that the 
effective nucleon mass becomes much smaller at high densities in the cases 
of PS9, PS10 and PS 11 than the other cases. This behavior is related to 
the small Mq in the parameters sets PS9, PS10 and PS11. From eq. (14), 
small Mq means large C v . In the parameter sets PS9, PS10 and PS11, the 
effects of the large C v overcome the effects of the negative Cq, and make the 
equations of state hard and M* small at high densities. It seems that the 
effects of the cu-nucleon couplings are more important than the effects of C 6 
at higher densities. However, it should be remarked that the higher terms 
of which can be neglected at the normal density may become important 
at higher densities, since $ becomes larger at that region. The effects of the 
higher terms at the higher densities are open questions. 



Figs. 6, 7 



4. Summary 

The results obtained in this paper are summarized as follows. 

(1) Using the effective Lagrangian which is constructed in the framework 
of renormalization group methods, we have studied the compressional prop- 
erties of the nuclear matter with the cutoff field theory, under the Hartree 
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approximation. 

(2) By including the corrections of the order (1/A') 2 , i.e., <£> 5 and $ 6 terms, 
into the effective Lagrangian, the physical results become to depend hardly 
on the value of the convenient cutoff A' which is conveniently introduced into 
the theory. 

(3) The compressional properties, the K-K', K-K c and K-K vs relations 
in the case of very large M$(= 0.9M) have the opposite features to the ones 
in the case of M * < 0.8M. 

(4) The coefficients of the effective Lagrangian were determined phe- 
nomenologically to reproduce the empirical data of K and K c . 

(5) In addition to the requirement of reproducing the empirical data of K 
and K c , if we require the stability of the solution of EOM of the a-meson at 
the limit $ — > oo, K = 200 ~ 300MeV are favorable. However, if the effects 
of the higher terms of $ may become important at large $ region, it is out 
of the scope of this approximation to discuss the stability at $ — > oo. If this 
condition of stability is not absolute in this approximation, the other choices 
of K are also possible. 

(6) K = 250 ~ 350MeV is favorable to reproduce K, K c and K vs simul- 
taneously. Naturally, this result is consistent with the result of the general 
discussion which does not depend on the details of the cr-meson potential 
[22]. 

(7) It seems that the effect of the cu-nucleon coupling is more important 
than the effect of C 6 at higher densities. However, the higher terms which 
can be neglected at the normal density may become important at higher 
densities, since $ becomes larger at that region. The effects of the higher 
terms at the higher densities are open questions. 

It is very interesting to reexamine the calculations of the meson self-energy 
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in the random phase approximation by including the (1/A') 2 corrections with 
the parameter sets obtained by this paper. It is now under the studies. 
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Table and Figure Captions 



Table I The sets of the empirical values of K, K c and K vs from the Table 
3 in ref. [19]. According to the conclusion in ref. [19], only the data in the 
cases of K — 150 ~ 350MeV are shown. ) 

Table II The parameter sets which reproduce the empirical value of K and 
K c in table I. K, K c , K 1 ', L, K sym and K vs are shown in MeV. (a) The sets 
with positive C 6 which reproduce the parameter sets of the values of K and 
the corresponding upper limits, mean values and lower limits of K c in table 
I. (b) The sets with negative C 6 which reproduce the values of K and the 
corresponding mean values of K c in table I. 

Fig. 1 The A'-dependence of K' . The solid, dotted and dashed lines corre- 
spond to the case of K =150, 250 and 350MeV, respectively, (a) The results 
with M* = 0.9 M. (b) The results with M* = 0.75M. (c) The results with 
M* = 0.6M. 

Fig. 2 K-K' relation in the case with A' =1.5GeV. The solid, dotted, dashed, 
dashed-dotted and bald solid lines correspond to the results with Mq = 0.5M, 
0.6M, 0.7M, 0.8M and 0.9M, respectively. 

Fig. 3 K-K c relation in the case with A' =1.5GeV. Each line has the same 
correspondence to Mq as in fig. 2. The crosses with error bars are the 
empirical data in table I, while the small solid triangles are the empirical 
data from Table IV in ref. [20] 
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Fig. 4 e' — de/d<& is shown as the function of $. The zero-point is the 
solutions of EOM of the a-meson. The solid, dotted, dashed, dashed-dotted 
and bald solid lines correspond to the results with the parameter sets, PS6, 
PS7, PS8, PS9 and PS10, respectively. 

Fig. 5 K-K vs relation in the case with A' =1.5GeV. Each line has the same 
correspondence to Mq as in fig. 2. The open circle at K — 250MeV, the open 
inverse triangles at K — 250 and 300MeV, the open triangles at K — 200 
and 250MeV, the solid circles at K = 150, 200, 250 and 300MeV, the solid 
square at K = 250MeV, and the solid inverse triangle at K = 350MeV are 
the results with the parameter sets PS3, PS2, PS5, PS1, PS4, PS6, PS7, 
PS8, PS10, PS9 and PS11, respectively. The crosses with error bars and the 
small solid triangles have the same correspondence to the empirical data as 
in fig. 3. 

Fig. 6 The binding energy Ej, is shown as the function of the baryon number 
density p. The solid, dotted, dashed, dashed-dotted, bald solid, bald dotted, 
bald dashed lines and the small solid circles correspond to the results with the 
parameter sets, PS1, PS2, PS3, PS4, PS5, PS9, PS10 and PS11, respectively. 

Fig. 7 The effective mass of nucleon M* is shown as the function of the 
baryon number density p. Each line and the small solid circles have the same 
correspondence to the parameter sets in fig. 6. 
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Set 1 


Set 2 


Set 3 


Set 4 


Set 5 


K 


150.0 


200.0 


250.0 


300.0 


350.0 


K c 


5.861 ±2.06 


2.577 ±2.06 


-0.7065 ± 2.06 


-3.990 ±2.06 


-7.274 ±2.06 


K 


66.83 ±101 


-46.94 ± 101 


-160.7 ±101 


-274.5 ± 101 


-388.3 ±101 



Table I 
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-76.22 



Table 11(a) 
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Table 11(b) 
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